Equilibrium correlations in charged fluids 
coupled to the radiation fleld 

Sami El Boustani, Pascal R. Buenzli,^ and Philippe A. Martin 
Institute of Theoretical Physics 
Swiss Federal Institute of Technology Lausanne 
CH-1015, Lausanne EPFL, Switzerland 



Abstract 

We provide an exact microscopic statistical treatment of particle 
and field correlations in a system of quantum charges in equilibrium 
with a classical radiation field. Using the Feynman-Kac-Ito represen- 
tation of the Gibbs weight, the system of particles is mapped onto a 
collection of random charged wires. The field degrees of freedom can 
be integrated out, providing an effective pairwise magnetic potential. 
We then calculate the contribution of the transverse field coupling to 
the large-distance particle correlations. The asymptotics of the field 
correlations in the plasma are also exactly determined. 

PACS numbers : 05.30.-d, 05. 40. -a, ll.lO.Wx 

1 Introduction 

Thermal states of non relativistic particles interacting by the sole 
Coulomb potential are known to provide an adequate description of 
many states of matter. The introduction of magnetic interactions be- 
tween the particles poses a novel problem since they are mediated by 
the coupling to the transverse part of the electromagnetic field. This 
immediately leads to consider the full system of matter in equilib- 
rium with radiation : the relevant theory becomes then the thermal 
quantum electrodynamics (thermal QED). 

In order to go beyond pure electrostatics without facing the full 
QED, a number of studies rely on the Darwin approximation. Dar- 
win has shown [1], [2] that one can eliminate the transverse degrees 
of freedom of the field within the Lagrangian formalism up to order 
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(c is the speed of light). A nice review of the derivation of the 
Darwin Lagrangian and a lucid discussion of its consequences can be 
found in [3] . The resulting Darwin Hamiltonian can be used to inves- 
tigate the equilibrium properties of the so called weakly relativistic 
plasmas; see the recent works of Appel and Alastuey [4], [5], [6] and 
earlier references therein. These authors have done a careful analysis 
of the domain of validity of the Darwin approximation and shown in 
particular that the predictions of the Darwin Hamiltonian on the tail 
of particle correlations in thermal states cannot be correct. Indeed the 
well-known Bohr-van Leeuwen theorem [7] asserts that classical (non- 
quantum) matter completely decouples from the radiation field. Thus 
the Darwin Hamiltonian, which treats the particles classically, should 
not predict any effect of the transverse field when used for thermal 
equilibrium computations. The Darwin approximation is, however, 
not deprived of any meaning in statistical physics. Indeed, the authors 
show in [6] that Darwin predictions about current correlations coin- 
cide with those of thermal QED in the restricted window of distances 
^part r ^ Aph, where Apart = h^J[3/rn is the de Broglie thermal 
wavelength of the particles and Aph = (5hc the thermal wavelength of 
the photons. But to determine the tail r S> Aph of the correlations in 
the presence of the radiation field, matter has to be treated quantum 
mechanically to avoid the conclusion of the Bohr-van Leeuwen theo- 
rem. The situation is similar to orbital diamagnetism in equilibrium, 
which is of quantum-mechanical origin. 

In this work, we consider equilibrium states of non-relativistic spin- 
less quantum charges coupled with the radiation field in the standard 
way (section 2). We shall, however, treat the field classically on the 
ground that the large distances r » Aph are controlled by the small 
wave numbers A; ~ - <C t^, implying fihuj-^ ~ ^ ^ 1. Hence only 
long-wavelength photons will contribute to the asymptotics which is 
expected to be adequately described by classical fields. The full QED 
model with quantized electromagnetic field will be studied in a sub- 
sequent work (see also comments in the concluding remarks, section 

8). _ _ 

Our main tool will be the Feynman-Kac-Ito path integral repre- 
sentation of the degrees of freedom of the charges. The Feynman-Kac 
integral representation has been widely used to derive various prop- 
erties of quantum Coulomb systems, in particular to determine the 
exact large-distance behaviour of the correlations; see [8], [9], and 
[10], [11] for reviews. In this representation quantum charges become 
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fluctuating charged loops (closed Brownian paths), formally analogous 
to classical fluctuating wires carrying multipolcs of all orders. These 
fluctuations are responsible for the lack of exponential screening in 
the quantum plasma and for an algebraic tail ~ of the particle 
correlations [12]. 

Adding an external magnetic field produces a phase factor in the 
Feynman-Kac-Ito formula, whose argument is the flux of the mag- 
netic field across the random loop. Correlations in the case of an 
homogeneous external magnetic field have been studied in [13] . When 
the particles arc thcrmalized with the field, the latter becomes itself 
random and distributed according to the thermal weight of the free 
radiation. The system can be viewed as a classical-like system of ran- 
dom loops immersed in a random electromagnetic field. At this point, 
the field degrees of freedom can be exactly integrated out by means 
of a simple Gaussian integral since the Hamiltonian of free radiation 
is quadratic in the field amplitudes. One is then left with an effec- 
tive pairwise current-current interaction between the loops which has 
a form similar to the magnetostatic energy between a pair of classical 
currents. For the sake of illustrating the basic mechanisms in a simple 
setting, this program is carried out in section 3 with particles obeying 
Maxwell-Boltzmann statistics. Appropriate modifications needed to 
take into account the particle statistics (Bose or Fermi) are given in 
section 7. 

In section 4 we apply the formalism to the determination of the 
asymptotic form of the correlation between two quantum particles em- 
bedded in a classical plasma. This simple model already illustrates the 
main features occurring in the general system. The effective magnetic 
interaction contributes to the tail, but its ratio to the Coulombic 
contribution is of the order of the square of the relativistic parameter 

(/?mc2)-l = (Apart/Aph)^. 

In section 5 we consider the generalization of the results obtained 
for two particles to the full system of quantum charges. The analysis 
relies on the technique of quantum Mayer graphs previously developed 
for Coulomb systems, and we merely point out the few changes that 
are needed to include the effective magnetic interactions. 

Field fluctuations in plasmas have been studied for a long time at 
macroscopic scales, much larger than interparticle distances; see [2], 
[14] and references cited therein. In section 6, we reexamine this ques- 
tion from a microscopic viewpoint and show that electromagnetic field 
correlations are always long ranged due to the quantum nature of the 
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particles. This is in disagreement with the prediction of macroscopic 
theories. Wc come back to this point in the conchiding remarks (sec- 
tion 8). However, in the classical limit, we recover the fact already 
observed in [14] that the long-range behaviour of the longitudinal and 
transverse parts of the electric field correlations compensate exactly. 

In section 7, we generalise the formalism developed in section 3 to 
include Bosc and Fermi particle statistics. This is done as usual by 
decomposing the permutation group into cycles and grouping particles 
belonging to a cycle into an extended Brownian loop. When this is 
combined with the Feynman-Kac-Ito path integral representation of 
the particles, the system takes again a classical-like form: a collection 
of Brownian loops immersed in a classical random electromagnetic 
field. At this point the physical quantities can again be analyzed 
in terms of Mayer graphs comprising pairwise Coulomb and effective 
magnetic interactions, as in section 5. 

The methods presented in this paper have been applied to the 
study of the semi-classical Casimir effect [15], [16]. 

2 The model 

We first consider the QED model for non-relativistic quantum charges 
(electrons, nuclei, ions) with masses Tn^ and charges contained in 
a box A G of linear size L and appropriate statistics. The index 7 
labels the S different species and runs from 1 to 5. The particles are 
in equilibrium with the radiation field at temperature T. The field 
is itself enclosed into a large box K with sides of length R, R L. 
The Hamiltonian of the total finite volume system reads, in Gaussian 
units, 

i=l '^^ i<j ' * i=l 

(1) 

The sums run on all particles with position r,, momentum pj, and 
species index 7^; Kvaiis(7j) ^j) is a steep external potential that con- 
fines a particle in A. It can eventually be taken infinitely steep at 
the wall's position, implying Dirichlet boundary conditions — i.e., van- 
ishing of the particle wave functions at the boundaries of A. The 
electromagnetic field is written in the Coulomb (or transverse) gauge 
so that the vector potential A(r) is divergence free and H^'^'^ is the 
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Hamiltonian of the free radiation field. The Coulomb gauge is usu- 
ally preferred for simplicity in situations where the particles are non- 
relativistic and high-energy processes are neglected [17]. It has the 
advantage to clearly disentangle electrostatic and magnetic couplings 
in the Hamiltonian. 

We impose periodic boundary conditions on the faces of the large 
box K ^. Hence expanding A(r) and the free photon energy in 
the plane-wave modes k = ^) gives 

A(r) = Y: 5('^);^(«kAe-^'^-^ + ay^^e^^n (2) 

kA 

i75^'i = ^^kaLakA (3) 

kA 

where aj^^^ and a^A a-re the creation and annihilation operators for 
photons of modes (kA), e^;^ (A = 1, 2) are two unit polarization vectors 
orthogonal to k, and Wk = ck, k = |k|. In (2), fl'(k), ^(0) = 1, is a 
real spherically symmetric smooth form factor needed to take care of 
the ultraviolet divergencies. It is supposed to decay rapidly beyond 
the characteristic wave number kc = mc/h (see [17], chap. 3). Since 
we are interested in the large-distance r — > oo asymptotics, related to 
the small-k behaviour k ^ the final result will be independent of 
this cut-off function. 

The total partition function 

Zl,r = TV e-^^^.« (4) 

is obtained by carrying the trace Tr = TrmatTrrad of the total Gibbs 
weight over particles' and the field's degrees of freedom : namely, on 

the particle wave functions with appropriate quantum statistics and 
on the Fock states of the photons. The average values of observables 
(Omat) = ^L^R^ (e~^'^^'^Omat) Concerning only the particle degrees 
of freedom can be computed from the reduced thermal weight 

PL,R yrad 

^Periodic conditions are convenient here. We could as well choose metallic boundary 
conditions. Since the field region K will be extended over all space right away, the choice 
of conditions on the boundaries of K are expected to make no differences for the particles 
confined in A. 



(5) 
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where Z^^^ = Tr^ad exp {—PHq^'^) is the partition function of the free 
radiation field, as follows from the obvious identity 

v 'I^mat {OmaX PL,r) 

{<^mat) - 7^ • iOj 

J-imat PL,R 

We shall perform the thermodynamic limit in two stages by first let- 
ting i? ^ CO. Then pi = liniR^oo PL,R defines the effective statistical 
weight of the particles in A immersed in an infinitely extended ther- 
malized radiation field. 

As discussed in the Introduction, in this paper we treat the elec- 
tromagnetic field classically. This amounts to replacing the photon 
creation and annihilation operators in (2) and (3) by complex ampli- 
tudes a^^^ and ajcA- In this case, the free field distribution factorizes 
out as exp {—PHr^l) = exp {—(3H}^) exp {—f3Hx), where 

Hr,l = Ha + H^"^, Ha = J2 (p^~ "^^^^^^) + C/pot(ri, 71, • • • , r„, 7n), 

i=l ^"^T' 



(7) 



and Upot is the total potential energy. Since the free radiation weight 
exp (— /Jifg'^'^) is Gaussian, A(r) = A(r, {okA}) can be viewed as a 
realization of a Gaussian random field, and the term Ha = -^^A({akA}) 
becomes the energy of the particles in a given realization of the vector 
potential having Fourier amplitudes {a^;^}. 
The partial trace (5) becomes, explicitly, 

where for a general function F{{ak\}) of the mode amplitudes {F)^^ 
denotes the normalized Gaussian average over all modes ^ 

(^)rad = n / ^ [pru.^e-^'^-\<^-^\'] F({akA}). (9) 

kA 

Note that the stability of Coulombic matter and the existence of 
thermodynamics for extended systems are assured if at least one of 
the species obeys Fermi statistics [18]. In the next section, merely 



■'The classical field is expanded as in (2) and (3) with dimcnsionlcss amplitudes ctkA- In 
fact there will be no h dependence arising from the field, as seen by changing everywhere 
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as a matter of simplifying the presentation, we compute the effective 
particle interactions defined by pi ignoring quantum statistics. In 
this case, Maxwell-Boltzmann statistics requires the presence of an 
additional short-range repulsive potential 14r(7i, Ir* — fjl) in the 
Hamiltonian (1) to prevent the collapse of opposite charges and guar- 
antee thermodynamical stability. The generalization to Fermi and 
Bose statistics will be given in section 7. 



3 The gas of charged loops and the ef- 
fective magnetic interaction 

We now introduce the Feynman-Kac-Ito path integral representation 
of the configurational matrix element (ri, r„|e~^^^|ri, r„) for 
the particles interacting with a fixed realization of the field. For a 
single particle of mass m and charge e in a scalar potential V'^^^{r) 
and vector potential A(r), we first recall that this matrix element 
reads [19], [20], [21] 



(rj exp 



X exp 



-/3 



(p-fA(r))^ 
2m 



/d^(s)-A(r + A^(s)) 

JO 

(10) 



f\s F"^*(r + A^(s)) - 
Jo 



VP' 



Here^(s), < s < 1, ^(0) = ^(1) = 0, is a closed dimensionless Brow- 
nian path and D(^) is the corresponding conditional Wiener measure 
normalized to 1. It is Gaussian, formally written as 



exp 



( 2/0 



ds 



ds 



')d[^(-)], 



with zero mean and covariance 



m)e{si)e{s2) = 5'^"(min(si, 32) - S1S2) 



(11) 



where ^^{s) are the Cartesian coordinates of ^(s). In this representa- 
tion a quantum point charge looks like a classical charged closed loop 
denoted by .7-" = (r,^), located at r and with a random shape ^(s) 
having an extension given by the de Broglie length A = fi^f^Jm (the 
quantum fluctuation). The magnetic phase in (10) is a stochastic line 
integral: it is the flux of the magnetic field across the closed loop. The 
correct interpretation of this stochastic integral is given by the rule of 
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the middle point; namely, the integral on a small element of line x — x' 
is defined by 



£de-f(0 = (x-xO-f(^), x-x'^0 



(12) 



We shall stick to this rule when performing explicit calculations.^ Note 
the dimensionless relativistic factor (/3mc^)~^/^ in front of the vector 
potential term. 

This is readily generalized to a system of n interacting parti- 
cles : The weight in the space of n loops = (ri, 71,^1), . . . ,J^n = 
(rn)7n)Cn) coming from the path integral representation of 
(ri,...,r„|e-^^A|r^^...^r^) is exp(-/3L'"(jri, A)) where 

n 

?7(jri , A) = ^ e^.e^^. VdJ^i, J^j) 

- /d^.(^) • A(r, + A,^.^,.(.)) 

(13) 



The matrix element (ri,...,r„|e '^^^ |ri, r„) is obtained by inte- 
grating exp(-/3[/(J'^i, A)) over the random shapes ^1, of 
the loops, as in (10). In (13), 

is the Coulomb potential between two loops, and for the sake of 
brevity, we have omitted the non electromagnetic terms 

n n 

Kr(^^, ^,) + E ^walls(^*) (15) 
i<j i=l 

corresponding to the short-range regularization and to the confinement 
potential. The vector potential term can be written as 



*Other prescriptions are possible for the path integral to correctly represent the quan- 
tum mechanical Gibbs weight in presence of a magnetic field. The Ito rule may be used 
when f is divergence free [20]. 
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—i Jdx A(x) • J7'(x) in terms of current densities associated with the 
Brownian loops : 



i=l 



m^^c^ Jo 



(16) 



If one interprets the (ill-defined) derivative A-y.d^j(s)/ds = Vj(s) as 
the "velocity" of a particle of charge e-y. moving along the loop ^j(s), 
the quantity e^.Vj(s)(5(x — rj — A^.^j(s)) corresponds to a classical cur- 
rent density. This is just a formal analogy. In subsequent calculations 
of stochastic integrals arising from (16), we will always use the math- 
ematically well-defined rule of the middle point (12). Moreover, such 
"imaginary time" currents appearing in the Feynman-Kac-Ito repre- 
sentation are not the physical "real-time" current observables. Our 
definition (16) also includes the relativistic factor (/9m^.c^)~^/^. 

A remarkable fact is that the transverse part of the field enters in 
exp{—f3U{J-i, ...,J^n, A)) as a phase factor linear in A and its Fourier 
amplitudes (contrary to the Hamiltonian (1) written in operatorial 
form). Since the statistical weight e'^^o" (3) is a Gaussian function 
of these Fourier amplitudes, it makes it possible to perform explicitly 
the partial trace over the field degrees of freedom in (8) according to 
the following steps : 



exp 



y"dx A(x)-J'(x) 



- / rad 



exp 



Y{ exp [i«A«kA + ^^kA«kA)] 



. kA 



rad 



2R^ ^ 

kA 



exp 



2 J (27r)3 



(17) 



The first equality is obtained by introducing the mode expansion (2), 
yielding 

n^x = a (^) ^ J(k) • ekA, J(k) = |dx e-i'- J(x) 

(18) 

The second equality results from (8), (9) and the Gaussian integral 
jd!a g-fe|a|2+i(«*a+na*) ^ ^-ig-fc-i |„|2 ^ ^ ^ whcreas the infinite vol- 
ume limit -R — > GO and the polarization sum have been performed in 
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the last equality. We have denoted by G^'^{\c) the covariance of the 
free transverse field : 

(19) 

{6^^(k) is the transverse Kronecker symbol). In (17) and throughout 
the paper, summation on repeated vector components = 1, 2, 3 is 
understood. In the configuration space, the asymptotic behaviour of 
Qiii^ (x) is obtained by approximating (k) ~ 1 in the inverse Fourier 
transform of G'^^(k): 

(20) 

Decomposing the total current (16) into the individual loop currents 
we see that the effective weight (17) takes the form 



(exp i/^ydx A(x)- J'(x) ^ =f[e^pl-f^W^{i,i) 



J / rad j=i 



exp j -^^^e-y^eJ.Wmi^,j) j , 

(21) 



where for two loops i = Ti and j = J^j we have introduced the loop- 
loop effective magnetic potential 

e^,e^,W^{i,j) = jd^jdy (/(J^,, x))*G'^''(x - y)/(JF,., y) = (22) 



/-r^e'^-^'-*-'-^) /Wf(si)e^''-H«»(^i) /"W(s2) e-*-S«^('^)G''^(k). 

J K^T^r Jo Jo 



As a consequence of Gaussian integration, one recovers pairwise inter- 
actions (22) between loops. The product in (21) contains the magnetic 
self-energies of the loops. 

It is pleasing and convenient that after averaging over the field 
modes, the energy of the system of loops becomes an exact and explicit 
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sum of pair potentials (and self-energies)^ 



(23) 



It is interesting to ask for the status of the partial density matrix 
(5) compared to that generated by the Darwin Hamiltonian poarwin 
g-/3HDarwin qi, more generally, if pl,r can be cast in the form pl,r 
Q-3H^fi gQjj^g tractable Hamiltonian ffcfr({Pi) ^j}) depending on the 
canonical variables of the particles. The answer to this last question 
is very presumably negative. Indeed the magnetic interaction (22) 
is a two times functional of the Brownian loops; namely, it lacks the 
equal-time constraint occurring in the Coulomb potential (14) (see the 
discussion before (26) below) necessary to come back to a simple oper- 
ator form by using the Feynman-Kac-Ito formula backwards. This is a 
well-known common feature of interactions resulting from integrating 
out external degrees of freedom [19]. 

The long-distance asymptotics of Wjn{i,j) as \ri — fjj oo is 
determined by the small k behaviour in the integrand of (22). Noting 
that Jq d^{s) =0 for a closed loop (Ito's lemma), one has 



/Wr(^) e^-^-H^.!^) ^ iX^^ l\^^s) k . Us), k ^ 0, 
Jo Jo 



(24) 



and thus 



W„.{i,j) ~ (25) 



fd^^{s,m{si) ■ VrJ /dej(^2)(e,(.2) • Vr,) G^^V, - Tj) 

Jo 



as |rj — Tjl oo. Upon using the asymptotic form (20) of (j'^^(rj — r^), 
it is clear that for fixed loop shapes and the decay of WYnihj) is 
~ |rj — rjl^-^. It is of dipolar type modified by the constraint imposed 
by the transversality. 

The Coulombic part (14) of the loop-loop interaction still decays 
as and deserves the following remark. From the Feynman-Kac 



^We omit again in (23) the non-electromagnetic terms (15). 
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formula the potential (14) inherits the quantum-mechanical equal-time 
constraint; i.e., every clement of charge e^. A^.d^,j(si) of the first loop 
does not interact with every other element ej.Xj.d^j{s2) as would be 
the case in classical physics, but the interaction takes place only if 
si = S2- It is therefore of interest to split 

V,{i,j) = V,i,,{i,j) + W,ii,j), (26) 

where 

Veiec{i,j) = [ dsi [ dS2 -. , , ^ , \ e / M (^'^) 

Jo Jo l^i + K^iisi) - rj - X^.^j{s2)\ 

is a genuine classical electrostatic potential between two charged loops 
and 



We(i,j) = [dSi [dS2{5{si-S2)-l) 

Jo Jo 



+ K^iisi) - rj - Xjj^j{s2)\ 

(28) 



is the part of Vc{i,j) due to intrinsic quantum fluctuations {Wc{i,j) 
vanishes if h is set equal to zero). Because of the identities 

/ dsi {d{si -S2)-1)= [ ds2 {5{si - S2) - 1) = 0, (29) 
Jo Jo 

the large-distance behaviour of Wc originates again from the term 
bilinear in and in the multipolar expansion of the Coulomb po- 
tential in (28) 

W,{i,j) ~ [dsi [dS2 {S{SI-S2)-1) {X^Msi) ■ VrJ {\^ji^2) " Vr,) , 

Jo Jo r 

(30) 

It is dipolar and formally similar to that of two electrical dipoles of 
sizes e^jA^.^j and e^jX^^^j. 



4 Two quantum charges in a classical 
plasma 

In order to exhibit the effect of the magnetic potential on the particle 
correlations, we consider the simple model of two quantum charges 
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and ei) with corresponding loops Ta = (ra,^^) and J^i, = (rb,^^) 
immersed in a configuration lo of classical charges, following section 
VII of [12] or section IV. C of [11]. According to (26) one can decom- 
pose the total energy as U {Ta , .^b , w) = ea^hWiTa ,J^h) + Uc\{J^a,H, 
where W{!FaT^b) = WdJ^a^^b) + ^^{J^a^^b) is the sum of the elec- 
tric and magnetic quantum dipolar interactions and UdiJ'ai^bi^) is 
the purely classical Coulomb energy (27) of the two loops J^a and J^b 
together with that of the particles in the configuration u. The corre- 
lation p{J^a-,^b) between the loops is obtained by integrating out the 
coordinates u of the classical charges : 

p{:Fa,Tb) = ^ fdu e-^"(^-^'"'^) =e-^^<'^''^(-^-^'')pci(^a,^fe), 

"cl J A 

(31) 

where Hd is the partition function of the classical plasma and Pc\{^a,^b) 
is the correlation of the two loops embedded in the plasma interact- 
ing with genuine classical Coulomb forces. In the latter quantity, the 
classical theory of screening applies so that effective interaction be- 
tween the loops decay exponentially fast ^. Thus one can approximate 
Pcl{^a^^b) in (31) by p{Ta)p{^b) up to a term exponentially decaying 
as |ra — r(,| oo. Furthermore, integrating p{Ta-,J^b) on the loop 
shapes leads to the following expression for the positional correlation 
of the quantum charges 

p(r„,r,) = |d(^j|d(^,) e-^^»^^^(^-^^)p(J-„)p(J-(.) + 0(e-^l'-''-"-''l) = 
= p„p;,-/3e„eb jBiuJmb) W{J^a,n)p{Up{ib)+ 

+ \p'elel jBiuJ^i^b) W^2(^a,^6)p(Op(^6) + - + C?(e-^l--^l) 

(32) 

Since W{J^a7^b) ~ ~ ^b]"^ (see (25), (30)), the above expansion 
in powers of W generates algebraically decaying terms at large sepa- 
ration. It is known that in a homogeneous and isotropic phase, the 
electric dipole part does not contribute at linear order [12], [11]. 
The same is true for the magnetic part. To see this, it is convenient 

^The usual Debye theory of screening has been rigorously shown to be valid at least at 
sufficiently high temperature [22]. 
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to write the linear Wjn term of (32) as 



= , IS^^'^-^'"'-'"'^ ir(k)t6(k)G^'^(k). (33) 

The stochastic ^„-hne-integral is now included in the definition of the 
tensor 



(34) 



and likewise for t^{^S-)- Since both the measure D(^„) and p(^a) 
variant under a rotation of in an isotropic system, ta(k) transforms 
in a covariant manner under rotations of k. Thus it is necessarily of 
the form ta(k) = fc^/a(|k|), implying the vanishing of (33) because 
of the transversality of G'*''(k). One concludes that the slowest non- 
vanishing contribution comes from the W"^ term in (32) 

The temperature-dependent amplitude A{I3) = AcdP) + ^mm(/9) + 
AcraiP) involves in principle electric and magnetic contributions from 
Wc^ and Wjij^, as well as a cross contribution from 2Wc^m- These con- 
tributions can be calculated explicitly at lowest order in h (or equiva- 
lently in the high-temperature limit /3 ^ 0) . The electric contribution 
in this limit is known to be [12], [11] 

AcciP) - h^^^PaPb. (36) 
240 maTni, 

To compute the magnetic contribution in the same limit, we write the 
quadratic term 

^JmMuJmb)Pi^b) wji^.,H) = 

^ /^/^"'^'''^'''^■^""''"''^ (rr(ki,k2))*rr(ki,k2)G'^'^(ki)G^ 

(37) 
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in terms of the tensors 

»1 rl 



rr(ki,k2) = hiUpiUl de(^i)/ <(S2) e-^^'^^^■^^is,)^-^x.^.,■Us,) 
J Jo Jo 

(38) 

and r^'^'^(ki, k2), defined likewise. As usual the behaviour at large 

distances is controlled by that of the integrand of (37) at small wave 
numbers. Expanding (38) at lowest order in ki and k2 gives 

rr(kl,k2) ~ [ma)piL)[d^iiisi)[dCis2) i-iXakl-Usi)) {-iXa\^2 ■ Us2)) 

J Jo Jo 

= -xikik^, [mjpiu fd^iiisi) fdc{s2) ea{si)ea{s2) 

J Jo Jo 

(39) 

and likewise for T^'^(ki, k2). One sees that because of the factor A^A^, 
the overall contribution in (37) will have a H'^ factor so that at this 
order we can neglect the quantum fluctuation in the density setting 
PiCa) ~ Pa independent of Thus the stochastic integral to be 
calculated becomes (appendix A) 



(40) 



leading to 



rr(ki,k2) ~ -^iS^'^k^ . k2 - k^k-o, 

rr(ki,k2) ~ -^(<5-ki • k2 - k^kl). (41) 

When this is inserted into (37) and summation on vectorial indices 
are performed, one finds the expression 



1 



^2^ 



'(42) 



with A = ^^^^^^af^^f^ _ fpj^g fjpg^ term in the large brackets gives a 

288mambC* ° 
rapidly decaying contribution since it involves the Fourier transform 

of the form factor ^^(k). The algebraic large-distance contribution 
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comes from the second term which reads, after Fourier transformation 
(approximating g{ii) ~ 1, k — ^ 0), 



A d^d,- . d^d,- 1 = A ^. (43) 

Finally one checks that there is no cross Coulomb-magnetic contribu- 
tion AcmiP) at the dominant order as a consequence of transver- 
sality (appendix B). So adding (36) and (43) gives the final result 



240 niamij 



1 + 



{(5maC^){l3mbc'^) 



1 



(44) 



as |ra — Tfil oo and at lowest order in h. One sees from (14) and 
(22) that the order of magnitude of the ratio VFm/K is (/3mc^)~^. In 
an electrolyte at room temperature T = 300K, this ratio is found to 
be « 10~^^. The magnetic correction to the correlation decay (44) is 
negligible in this case. 

5 Particle correlations in the many- 
body system 

We apply the formalism developed in section 3 to the determination 
of the large-distance decay of the particle density correlations in the 
more general case where all particles are quantum-mechanical, but 
still obeying Maxwell-Boltzmann statistics. 

We show hereafter that the algebraic decay of the (truncated) 
particle density correlations 

pT{7a,ra,7b,n) \ — — k-, |r„-r6|^oo (45) 

found in the absence of the radiation field [9], [11] is not altered, but 
that the coefficient Aab{P,{Pj}) contains in addition small magnetic 
terms of the order (/9mc^)~^, as in (44). As an illustration, we give the 
lowest order contribution of this coefficient with respect to Planck's 
constant h. 

By the Feynman-Kac-Ito representation, the full system composed 
of quantum point charges coupled to the radiation field has reduced 
to a classical-like system of extended charged loops = (r,7,^) for 
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which all the methods of classical statistical mechanics apply. The 
only novelty comes from the additional magnetic potential Wm. In 
the following, we merely summarize the arguments since they are es- 
sentially the same as those found in [9], [11] when no radiation field 
is present. 

As usual, we express the truncated two- loop correlation p^[Ta,J'b) 
= p{J^a)p{,^h)h{J^ai^b) in terms of the loop Ursell function h{J^ai^b)- 
The latter function can be expanded in a formal diagrammatic Mayer 
series of powers of the loop densities p{T). One needs to resum the 
long-range part of the Coulomb potential V^, which is responsible for 
the non-intcgrability of the Mayer bonds 

f{Ti,J^j) = exp{-f3e^^e^.[Vc{J^i,J^j) + WmiJ^i,J^j)]) - 1 at infinity. Us- 
ing the decomposition (26) we resum the convolution chains built 
with the purely electrostatic long-range part Vc\cc{^, into a Debye- 
Hiickel-type screened potential ^eiec(-^! ^')- Then reorganizing the di- 
agrams leads to a representation of the loop Ursell function by terms 
of so-called prototype diagrams, built with the two kinds of bonds 

F(.F,.F') = -/3e-,ey$eiec(^,^'), (46) 

F^(^,^') = e-/?-7ey[*elec(^,^')+W^{^,^')] _ 1 +/3e^ey$elec(^,^'), 

(47) 

where we have defined W = + W^a as in section 4. ^ 

The potential $eiec(-^j -^O has been studied in [23]. It corresponds 
to the term n = of the full quantum analog of the Debye-Hiickel 
potential given by formula (89) of [23]. This contribution n = is 
shown to be decaying at infinity faster than any inverse power of |r— r'| 
(see formula (58) of [23], and the comment following it). 

The asymptotic decay of the two-particle correlation pTila, Ta; 7fe; fft) 
is inferred from that of the loop correlation p-Y{ J~a; ^b) by integrating it 
over the Brownian shapes and ^(,. The bond F is rapidly decreasing, 
and the asymptotic decay of is dominated by the dipolar decays 
of Wc and Wra : F^{J^,P) ~ -(3e^e^,W{T,P) as |r - r'| ^ oo. We 
further extract this dipolar part from F^ and define the bond 

FR(jF,jr') = F^{J^,T')+l3e^e^W{T,T') 

~ \\J3e-ye^>W{J^,J^')f = 0(|r - r'p^) (48) 

^Strictly speaking, the short-range repulsive potential needed in the framework of 
Maxwell-Boltzmann statistics would arise here in the exponent of (47). It has no im- 
plication in this discussion about long-range behaviours, and we simply omit it. 
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and work now with the three bonds F, F^, and W. ^ 

To find out the slowest-decaying diagrams, we write the truncated 
two-loop correlation pT{^a, ^h) in an exact Dyson series of convolution 
chains involving W and H : 

+ p'^{Ki.Wi.Ki.Wi.K){J^a,H) + - (49) 

where H denotes the sum of the diagrams that remain connected un- 
der removal of one VF-bond and K{J^i,J^2) = p{,^i)p{^2)H{J^i,J^2) + 
5{Ti,T2)p{^i)- This topological constraint ensures that H decays at 
least as r~^. The series (49) is conveniently analysed in Fourier rep- 
resentation with respect to — rfe. After expanding W into the sum 
Wc + Wra, we have three types of chains : pure Wc or chains 
and mixed VFc, Wm chains. It is shown in [9], [11] that the contribu- 
tion of pure Wc chains to the particle correlation /9t(7oj ^a, 7;,, rj) = 
/D(^a) /D(^b) pT(-^a,-^6) decays strictly faster than o(|ra - rfel"*^). 
^ We show below that all other chains containing W^a bonds vanish 
identically as the consequence of transversality. This implies that the 
longest-range part of the correlations originates from the function H 
in the first term of the right-hand side of (49), hence the result (45). 

A chain mixing Wc and W^ bonds must have at least one element 
Wc -k K -k Wra 01 W^ * K k Wc- In Fourier space, one can write, from 
(28) and (22), 

{WckKkW^){ja,L,lbM,^) = (^(^« - ^i) - l)^e^'^-^-''^''(^'') 

X [T''^{k,si)G'"'^''^(k)] [^dC'^'{sb)e-'^-^^"'^^^'^\ (50) 

Jo 

where 

T'^Hk,.i)=^/"D(^i);^/"D(6)e-^'^-^-i^i(^i)if(7i,^i,72,6,k) 

71 72 

X /de2'(s2)e"'-^^2^^(^2) (51) 

**In [11], the bond F is further decomposed into a multipole expansion. Our bonds F^ 
and F^ differ formally from their bonds F\ and F\ only by the inclusion of the magnetic 
contribution VFm into W . 

^In this proof, only the invariance of H under rotations is used, which also holds when 
the magnetic potential is included. 
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and if (71,^1, 72,^2 5 k) is the Fourier transform of K{Ti,T2) witli re- 
spect to ri — T2. As the measures D(^i) and D(^2) and the func- 
tion -f^ (71,^1, 72,^2) k) are invariant under spatial rotations, T'^'^^k, si) 
transforms as a tensor, implying that it is necessarily of the form 
T^'^i\<i,si) = k'^^ a{k,si) for some rotationally invariant function a of 
k. Using fc^C^(k) = one deduces immediately that (50) vanishes. 
The case of -k K -k Wc is similar. To see that there are no chains 
containing only Wm bonds in pTi'Ja, ^a, lb-, ^h)-, it is sufficient to notice 
that the integrated root element J'D(^o) K ★ also involves a fac- 
tor [r^'^{\i)G'^^'^''{k)] (for another function T^'^{k) transforming in a 
covariant manner), and thereby vanishes for the same reason. 

The graphs that do contribute to the coefficient Aab{fi , {p-^}) of 
(45) are those of H that contain bonds with algebraic decay : namely, 
and W . To select the lowest contribution in /i, one notes first 
that W is at least of order fv^, as seen in (25), (30) which correspond 
to the lowest-order terms in the multipolar expansions of and 
Wra.. (Higher-order multipoles generate higher powers of the de Broglie 
wavelengths.) Since $eiec is rapidly decreasing, the algebraic part of 
is of order and is given by \[[3e-^e^iW{T ,J^')Y , as in (48). 
Thus, up to order fi^^ graphs with an algebraic decay can contain 
only one bond W , two bonds W , or one bond F^ belonging to paths 
connecting the two root points. If there is a single such link W, by 
the topological structure of H there exists another path connecting 
the root points made of the more rapidly decreasing bonds F and 
F^. Hence the whole graph has a decay faster than . If there are 
two W bonds in between the root points, as each of them is of order 

all the other bonds and vertices can be evaluated in the classical 
limit h ^ Consequently, at least one of the extremities of either 
bond W is attached to a purely classical part of the graph, which is 
independent of the Brownian shapes. We call such a point a classical 
end of W . At such points, integration over the Brownian shape of 
the loop "kills" the decay of W (see Appendix C), leading to an 
overall decay faster than r^^ . Finally, at order /i^, the only graphs 
that contribute to (45) arc constituted by a single F^ bond linked to 
the root points by purely classical subgraphs. The sum of such graphs 
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contributes to the particle correlation in the large-distance limit as 



7l>72 



(52) 



D(^2 



/3e^i e^2 VF'^'P (71 , ^ 1 , 72 , 6 , i-a - r^) 



where W^'^ = W/'p + W^'^'p is the sum of the dipolar parts (30) 
and (25) of Wc and W^, and ra^(7a,7i,r) is the classical truncated 
two-point density correlation (including coincident points). The func- 
tional integrals in (52) have been calculated in section 4, see (37)-(44), 
yielding the final result 

prha, ra, lb, n) ^^^^ X] Z'^*' "'T(7a, 7i, r) / dr 71^(72, 76, r) 



7i!72 

^11 ^11 



1 + 



1 



^6 



(53) 



as |ra — r^l ^ 00 and at lowest order in h. To this order, the only 
difference with (44) is the occurrence of the classical correlation func- 
tions n^, a manifestation of the fact that in the quantum many-body 
problem, every pair of particles contribute to the tail of the correla- 
tion function. This generalizes the result of [12], formula (5.12), to 
the inclusion of the magnetic interactions. 

As a final comment, we observe that the inclusion of the transverse 
degrees of freedom of the field does not modify the charge sum rule in 
the system of loops and hence it also holds for the charge correlations 
in the particle system. This sum rule reads 



dr/D(e; 



y-^ e^PT(-^,-^i) 



7 



= — e- 



■71- 



(54) 



It states that the charge of the cloud of loops induced around a fixed 
loop J^i exactly compensates that of Ti. The proof can be carried 
out word by word as in [23], section 6.1.2. It relies exclusively on the 
long-range part r^^ of the Coulomb potential Vc and is not altered by 
the presence of the magnetic potential Wra- 

6 Transverse field correlations 



A characteristic feature of charged systems is that longitudinal field 
correlations always remain long ranged in spite of the screening mech- 
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anisms that reduce the range of the particle correlations. It has been 
established on a microscopic basis that the correlations of the longi- 
tudinal electric field Ei behave as [24], [25] 



1 



x-y| 



27r 
' 3 



j dr |r|^5(r) 



|x - y| ^ GO, 
(55) 



where S{r) is the (classical or quantum-mechanical) charge-charge cor- 
relation function. 

In order to obtain the correlations of the transverse fields we first 
consider correlations {A^^{yi)A'^ {y))^ of the vector potential at free 
points X and y in space. These correlations are easily obtained by 
functional differentiation, adding to the original Hamiltonian (1) a 
coupling to an external current Sext{^) 



ff.,H(J„.) = ff.,H-^/dxJ„.(x).A(x), 

SO that 



(56) 



1 
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■ InTr Q~^^L,R{Je:i!t) 



3 ext— 

(57) 



Decomposing as in (7) one can write 
(^'^(x)^''(y))T = 



In TVmat ( e 



-PHa JP f dx J'ext(x)-A(x) 



rad 



ext —0 

(58) 



Using the Feynman-Kac formula as in section 3 one sees that the only 
modification in (17) is the replacement of the loop current »7(x) by 
the total current^*^ 



Jtot(x)=J'(x)+J'ext(x). 



(59) 



a consequence of the imaginary coupling constant in the Hamiltonian (56), the 
total current is real, so that we can still apply the Gaussian integration formula used in 
(17). 
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The Gaussian integration on the field variables replaces (17) by 

X [{jn* r + (JTxt)* + (^'^)* ^ext + U^t)* ^ext] (k)}- (60) 

Therefore, from (60), functional differentiation with respect to ,7ext 
according to (58) produces two terms 

(^'^(x)^'^(y))T = {A^{yi)A\y))\ + (A'^(x)A'^(y))?-^ (61) 

The first contribution (written in Fourier form) 

-^rV"^^^)^ r^oo, r = x-y, (62) 

arises from the part quadratic in Jext in (60) . It describes the thermal 
fluctuations of the free field, and in view of (20), decays as r~^. The 
second term, coming from the part linear in ,7ext , 

{A^{:x.)A''{y))T' = 

= ~ /(|^e*-y^e^^'-''G'^'^(k)G-(k') (^^(k)J-(k'))T, (63) 

represents the modification to the free-field fluctuations caused by 
the presence of matter. It involves the loop current correlation func- 
tion (J''^(k) J"^(k'))T where the average is taken with respect to the 
thermal weight (23) for the statistical-mechanical system of loops. 
Expressing the currents »7"(k) = f dJ^ j(jr, k)/5(.F) in terms of the 
density of loops p(J^) = Y^i^{^-,^i) (see (16)), one can write this 
current correlation in terms of the loop density correlation function 
^T(7i,Ci,725^25k) (including the contribution of coincident points) : 

(J"(k)J^(k'))T = (27r)^<5(k + k') 

X E hi^i)hi^2) T'^(7i,^i,k)(r-(72,6,k))V(7i,6,72,6,k). 

(64) 

The (5(k -|- k') is the manifestation of the translational invariance of 
the state, and we have set 

r-(7,,C,,k) = ^il= /dC(sO e^^-.'^-^*^^*). (65) 

y pm^^c Jo 
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When (64) is introduced into (63), one obtains the final form 

(A'^(x)A^(y))?J^t = - e*i'-("-y)G'^'^(k)G'^-(k)Q'^^(k), (66) 

where Q'^'^ (k) is the tensor 

Q--{k)=J2 y"D(ei)y"D(6)r-(7i,6,k)(r-(72,6,k))V(7i, 6,72,6, 



71,72 



(67) 



To obtain the long-distance behaviour of this correlation we examine 
the integrand in (67) at small k. Because of isotropy, the tensor 
Q'^'^(k) transforms covariantly under the rotations, and thus is of the 
form 

Q^^(k) = a{k)6^^ + h{k)k''k^, k = |k| (68) 

The term k^k'^ does not contribute to (67) since G^'^(k) is transversal. 
Because of Ito's lemma, '?''^(7i,6, k) is linear in k as k — > 0, implying 
a{k) = a + o{k)]. Hence, using 5f;'(k)5J7(k) = 5f;'(k) one finds 

G'^^(k)G'^^(k)Q^^(k) = 47ra ^ C(k)[l + o{k)] = 47r a G^'^ (k) [1 + o{k)] 

(69) 

as fc ^ 0. This shows that (A''(x)A''(y))J['^' has the same type of 
decay as the free field part (62) with a modified amplitude. Summing 
up the two contributions (61) gives 

(A^(x)yl-(y))T - ^ [s^"" + ^) Q - 47ra) , r ^ oo. (70) 

For B(x) = V X A(x), one finds, from (70), 

(5''(x)5''(y))T ~ (^M^.^) - 47ra) , r ^ oo. (71) 

The constant a = a{h, j3, p) embodies the effects of matter on the 
transverse field fluctuations. It has a relativistic factor (mc^)~^ and 
vanishes in the classical limit /i ^ (in accordance to the Bohr— van 
Leeuwen decoupling) as 0(^'^) (see Appendix D). 
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Et(x) 



In order to find the correlations of the transverse electric field 
1 dA{x,t) 



dt 



t=o 

1/2 



kA 

(72) 

we couple the latter to an external polarisation Pext (x) , 

HL,R{Ve.t) = Hl,R - ^/dx Pext(x) • Et(x), (73) 

and proceed as after (56). This amounts to replacing everywhere 
J7'cxt(k) by ikPextO^) so that the right-hand side of equation (60) is 
changed into 

X [(J")* J" - ik (V^S r + it (3")' P« + (P£:t)' P«] (k)}. 



As Pcxt(r) and »7(r) are real, 'P*xt(k) = ■Pext(— k) and likewise for J . 
From the change of variable k i-^ — k, one sees that the second term 
in the integrand in (74) is exactly compensated by the third term. 
Only the term quadratic in T'ext remains, which is responsible upon 
functional differentiation for the thermal fiuctuations of the free field, 
as in (62). Hence, the correlations of the transverse part of the electric 
field are decoupled from matter and one finds 

(x)£;ay))T = (£;rWK(y))^ - (^m^^^) \. ^ ^ oo. (75) 

The asymptotic correlation of the complete electric field E(x) = Ei(x)+ 
Et(x) follows from (55) and (75) (one can check by similar calculations 
that the cross correlation (£^,'*(x)£'^(y))T vanishes identically) : 

(i?^(x)E^(y))T = (i?f(x)i?r (y))T + (£;nx)sny)>T 



oo. 
(76) 



In the classical limit, S(r) satisfies the second-moment Stillinger- 
Lovett sum rule [25] — ^ Jdr |rpS'(r) = 1/13. Hence, the asymptotic 
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longitudinal electric field correlations in the matter are exactly com- 
pensated by those of the free radiation part, as noted in [14]. However, 
this no longer holds for quantum plasmas. As an illustration, for the 
quantum one-component plasma (jellium), one has [26] 

/dr M^S(r) = "f^coth (^) = i + 1(5^) % O(a'), 

(77) 

where u)p is the plasmon frequency. The long range of the electric 
field correlations is thus a purely quantum- mechanical effect. These 
findings are further discussed in the concluding remarks (section 8). 

7 Bose and Fermi statistics 

In this final section we introduce the needed modifications when the 
Fermionic or Bosonic particle statistics are taken into account. 

The Bose or Fermi statistics of the particles can be incorporated 
in the formalism following the same procedure as described in [8], [11] 
(section V). The matrix elements of (8), which is still an operator 
depending on the particle variables, are to be evaluated with prop- 
erly symmetrized (antisymmctrizcd) states. When combining the de- 
composition of the permutation into cycles with the Feynman-Kac-Ito 
path integral representation this leads to generalize the previous loops 
T = (r, 7,^) to Brownian loops C = (g, R, 7, X) that carry q particles 
(a set of particles labeled by indices belonging to a permutation cycle 
of length q) . The loop is located at R and has a random shape which 
is a Brownian bridge X(s), < s < q, X(0) = X((/) = with zero 
mean and covariance 

We merely give the final formulae since all steps are essentially iden- 
tical as those presented in the above mentioned works. 

The grand canonical partition function of the particle system, 
with the field degrees of freedom integrated out, takes the following 
classical-like form in the space of loops 

Sa = ^ -: HdCi z(A)exp ( - pU{£i, . . . ,Cn)). (79) 

n=0 i=l 







mm 1 — 




V Q 





£1 S2 



(78) 
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Integration on phase space means integration over space and summa- 
tion over all internal degrees of freedom of the loops : 



ldC--- = /dR^f; IdCX).... (80) 



7 9=1 



U {Ci, . . . , Cn) is the sum of the pair interactions between two different 
loops e^^e^.[VciCi,Cj) + y^mi^i,^j)] with 

r<}i r<}j 1 
Vc(A,>Cj)= / dsi/ ds2 S{si -h) , , . — ^ X ^ ( \\ 

Jo Jo \ti4 + ^'Ti-^iKSl) — IXj — A^jJ^j(S2)\ 

(81) 

the Coulomb potential, and 

X /^dXf(si)e*-H^*(^i) f\x^As2)e-'^'^^^^'^'^'^ G>'^{\^) 
Jo Jo 

the effective magnetic potential obtained after integrating out the field 
variables. Here s = s mod 1 and 5{s) = J2'^=-oo e^^^"'^ is the periodic 
Dirac function of period 1 that takes into account the equal time 
constraint imposed by the Feynman-Kac formula. Finally, the activity 
z{Ci) of a loop 



(27rg,A2J3/2 

(83) 



z(£,) = ™_ exp(-/?[Z^,eif(A) + Vwaiis(A)]), = e^"-' 



incorporates the chemical potential /x-y. of the particle, the effects of 
quantum statistics (rj^. = 1 for bosons and rj^. = —1 for fermions), 

and the internal interaction Useii{jO-i) = — 2^(Vc + yVm){^i,^i) of the 
particles belonging to the same loop (omitting the infinite Coulomb 
self-energies of the point particles). The addition of the magnetic po- 
tential Wm is the only modification compared to the formalism previ- 
ously developed for pure Coulombic interactions. Maxwell-Boltzmann 
statistics and the potentials (22) and (14) of section 3 are recovered 
when only single-particle loops {q = 1) are retained. 

At this point, due to the classical- like structure of the partition 
function (79), methods of classical statistical mechanics can be used 



26 



in the auxiliary system of loops, in particular the technique of Mayer 
graphs, as in section 5. The statistics of the particles affects the 
coefficients of the tails of the density and field correlations, but not 
their general forms (45), (71) and (76). 

8 Concluding remarks 

The Fcynman-Kac Ito path integral representation of the Gibbs weight 
enables one to integrate out the (classical) field variables. This yields 
an exact pairwise magnetic potential in the space of loops, which is 
asymptotically dipolar. It generates small (0((/?mc^)~^)) corrections 
to the tail of the particle correlation due to the pure Coulombic inter- 
action. 

A word is necessary about spin interactions that have not been 
included in the above discussion. Spin-1/2 electrons give rise to the 
additional term — ^^^iLif j • B(rj) in the Hamiltonian, with B(r) = 
VA A(r), V = gs the gyromagnetic factor, and a the Pauli matri- 
ces. Using spin coherent states [27], a functional integral representa- 
tion of the Gibbs weight can be constructed including the coupling of 
the spins to the field. Since this coupling is linear with respect to the 
vector potential, Gaussian integration on the field variables leads again 
to an effective spin-spin interaction VFs(^)j) and effective cross inter- 
actions Wsm(^,j) and WrDs{i-ij) between spin and orbital magnetism; 
details can be found in [28]. One finds that these interactions are of 
dipolar type ~ r~^, r ^ oo and they have to be added to the magnetic 
potential W^{i,j). In a homogeneous and isotropic phase, the spin 
interaction terms contribute to the tail of the particle correlations 
with the same amplitude °^ mlc^ numerical factors, as that 

found in section 4 in the case of the magnetic potential Wj^. 

Regarding the electric field correlations in the plasma, we also find 
that they have an algebraic decay of dipolar type. This is in dis- 
agreement with the macroscopic calculation presented by Landau and 
Lifshitz [2], §88, based on linear response theory and the fluctuation- 
dissipation theorem. Indeed, the electric field fiuctuations obtained in 
this theory are short ranged (exponentially fast decaying) : unlike in 
our calculation, the algebraic parts of the longitudinal and transverse 
correlations compensate exactly in the Landau and Lifshitz formu- 
lae [29]. Understanding the relation between our strictly microscopic 
approach and the macroscopic theory of field fluctuations is an open 
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problem. 

Let us, however, briefly point out some differences between the 
two approaches. The microscopic approach involves all length scales, 
whereas Landau and Lifshitz assume that the correlations of the po- 
larisation arc local (S correlated in space) and thus deal with a wave- 
number-independent dielectric function e{uj). Taking into account 
the wave-number dependence, it is likely that e(k, w) has terms non- 
analytic in k, reflecting the fact that Coulombic matter has alge- 
braically decaying correlations. In fact, for the jellium model, the 
static dielectric function e(k,u = 0) has a singular term ~ |k| in its 
small-k expansion [30]. It is possible that in a non-local generaliza- 
tion of the Landau-Lifshitz theory such singular terms also generate 
power-law decays of the field correlations. Furthermore, the magnetic 
permeability is usually set equal to that of the vacuum, thus ignoring 
the magnetization fluctuations, whereas in our calculation the latter 
are properly included. 

We stress again that the results of this paper hold when the elec- 
tromagnetic field is classical, which is supposed to correctly depict 
the small-wave-number regime, as said in the Introduction. Hence, 
the occurrence of the Planck constant originates exclusively from the 
quantum-mechanical nature of matter. If the field is quantized, we 
can, however, not exclude an interplay between /imattcr and ^gcldj 
which could lead to a modification of the asymptotic formulae pre- 
sented in the paper. 
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Appendix A 

To establish (40) according to the middle point prescription (12) one 
has to evaluate the rotationally covariant tensor 

[b{0 fd^is) fdCit) e{s)m = (84) 
J Jo Jo 

n,m 

where kn = ^ and Im = j^- Setting C{s,t) = 5'*''(min(s, t) — st) (see 
(11)), one has 

~ n m™oo 4 ~ ^ {^rn ^m~^l) ~ C {kn~\^ Irn) + C (^n, — ^, ~ ^)] 

X {kfi, Im) C (^kji, Im — ) -\- C (^kfi — ,/rn) + C(fen TJj^m m)^ ' 

A2 = ^Ihn^ \ {C {kn, kn) + C [kn, A^n" ^) - C (fc^ - 1, fc„) - C - 1, fc„ - 1)] 

X I^C {Im ) Im ) ~l~ C' (/jn , Im m^ ^ (^^in m ' ^ n ' n ) ] ' 

(85) 

This results from the application of Wick's theorem to the Gaus- 
sian average (84) with covariance (11). Expanding C [k^ — ^,lm) = 

C {kn,lm)-^{dlC) {kn,lm) andC{kn,lm - ^) = C {kn,lm)-^{d2C) (knjm) 

and taking the limits n, m — ^ 00 gives 

^1 = [ds [dt C{s,t){did2C)is,t) = 4, 
Jo Jo J-^ 

An = (ftC)(«,t) (dtC){s,t) = -i, (86) 

hence the result (40). 
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Appendix B 



From (22) and (28) the cross Coulomb-magnetic term is 




S2) - 1) 



(87) 



■1 



where 




-aoki-{„(si) 




,-jAak2-€a(s) 



(88) 



Because of the rotational invariance of D{^g)p{^^), averages of odd 
powers of ^„ vanish. This implies that in the small-ki , k2 expansion 
of Ha{ki,k2,si) only odd monomials in ki,k2 occur: 



where 03(ki,k2) represent monomials of order 3 in the components 
of ki,k2. The same holds for if^(ki, k2, S2). Since k^Gi''' {k2) = 
by transversality, one concludes that the term (87) decays at least as 
\ra-rb\'^. 

Appendix C 

If point i in Wc{i,j) or Wm{i,j) is a classical end, there is no other 
dependence at this point than that arising from these bonds. In the 
asymptotic formula (30) for Wc, this dependence is linear and thereby 
vanishes upon the space-inversion invariant D(^j) integration. In the 
case of Wm, from formula (22), the D(^j) integration yields the factor 



J Jo 

because of covariance under rotation. Hence, this contribution van- 
ishes as a consequence of transversality k'^G'^^{k) = 0. 




•^„(s))+03(ki,k2) 



= const X + 03(ki, k2), 



(89) 




(90) 
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Appendix D 



An explicit expression for the constant a = a(h, (5, p) follows from 
taking the trace in equation (69) and using (67) expanded for small 
k. This yields 



7>7 

"1 rl 



de{s) d^is') (k • ^{s)) (k • e'(s'))C(k) utM, 7',r, k = 0), 
Jo 

(91) 

where k = k/Zc. As A-yA^/ is of order h^, at lowest order in h one can 

set ^ = in the correlation function. The latter becomes independent 
of the quantum fluctuations ^, ^' and reduces to the density correlation 
function of the corresponding classical system. The remaining func- 
tional integrals, of the type /D(^) d^^^{s)^'^{s), vanish identically. 
The terms of order O(^) in nx arc necessarily linear in £ or They 
do not contribute to a since averages of odd powers ^ or ^' are zero, 
implying that there are no /I'^-terms in a. We thus conclude that a is 
0{h% 
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